Polarized gluon studies with charmonium and bottomonium at LHCb and AFTER 
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Recently it has been put forward that linearly polarized gluons inside unpolarized protons affect 
the transverse momentum distribution of final state particles in hadronic collisions. They lead to 
a characteristic modulation of the differential cross section in Higgs production and to azimuthal 
asymmetries in, for instance, heavy quark pair production. Here we study the effect on charmo- 
nium and bottomonium production in hadronic collisions, such as at LHCb and at the proposed 
fixed target experiment AFTER at LHC. We focus mainly on the scalar and pseudoscalar quarko- 
, nia, 77c, XcO, ?7t>, Xbo, which allow for an angular independent investigation. Within the framework of 

. transverse momentum dependent factorization in combination with the nonrelativistic QCD based 

color-singlet quarkonium model, we show for small transverse momentum (q^ <C 4Mq) that the 
^\ | scalar and pseudoscalar quarkonium production cross sections are modified in different ways by 

CN| . linearly polarized gluons, while their effects on the production of higher angular momentum quarko- 

nium states are strongly suppressed. Comparisons to Xc2 , Xb2 production can help to cancel out 
uncertainties. Together with the analogous study in Higgs production at LHC, quarkonium produc- 
tion can moreover be used to test the scale dependence of the linearly polarized gluon distribution 
over a large energy range. 

PACS numbers: 12.38.-t; 13.85.Ni; 13.88.+e 
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\Q ■ It has been recently pointed out [1, 2] that the study of bound states of heavy quarks (quarkonia) with even charge 
CO | conjugation, such as r\ c _b and Xc,b, produced in high-energy hadronic collisions, provides a useful tool to access the 
00 ■ unpolarized gluon distribution function. These quarkonia are produced by fusion of two gluons at scales of the order 
of their mass, thus large enough to justify the use of perturbative QCD. At leading order, they originate from a 2 — » 1 
partonic reaction, with no additional gluon emission in the final state. For this reason event rates should be high 
and, in analogy to the Drell-Yan process, they have very simple kinematics, with gluon momentum fractions directly 
related to the rapidity of the observed quarkonium state. This is in contrast to the production mechanisms of J/ip, 
which has earlier been suggested to pin down the gluon content of the proton [3, 4]. C = + quarkonia should be 
more reliable as gluon probes, because one expects neither large QCD corrections, nor the many open questions and 
theoretical uncertainties that affect the predictions of J/ip and T production rates [1, 2]. Sec Rcf. [5] for a recent 
review on this topic. QCD corrections to C = + quarkonium production have been calculated in Refs. [6, 7]. 

Hadrons resulting from 2 — > 1 scattering processes typically have a small transverse momentum, determined by the 
transverse momenta of the partons in the initial state. Hence, while being mostly lost down the beam pipe at collider 
facilities like the Large Hadron Collider (LHC), they could be detected by forward detectors (LHCb) or in fixed target 
experiments (AFTER@LHC). Here we do not aim to make a study of the optimal quarkonium decay channels for 
these experiments, but according to Ref. [8], the forward detectors and powerful particle identification at LHCb allow 
for a complete study of all charmonium states through their pp decays. C-even charmonium decaying to 4><j) provides 
another channel suitable for LHCb, because of the large branching ratio of cf> — > K + K~ and the clean signature of 
two narrow (f> signals [8]. Furthermore, at AFTER, a proposed fixed-target experiment which would utilize the 7-TeV 
proton beam of the LHC extracted by a bent crystal, yielding collisions at yfs ~ 115 GeV, one could study the r\ c ^ 
states in the 77 channel and the Xc,f>2 states through £+^"7 decays [1, 2]. For discussions of promising channels to 
study 77f, production, cf. e.g. Refs. [9-11]. 

The states that we will mainly focus on in this paper will be the even charge conjugation, scalar and pseudoscalar 
quarkonia, i.e. J PC — ±+ states, or in terms of spectroscopic notation ( 2S+1 Lj) the 1 Sq states r? c ,?76 and 3 Pq states 
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XcOiXbo- The XdiXbi states would require a different treatment, because they suffer from the same problem as other 
vector states, such as the J/i/j, namely that because of the Landau- Yang theorem, their production from two gluons 
requires either an additional gluon radiated off or inclusion of the off shellness of gluons. Although off-shell gluons 
contribute significantly to Xd production according to Rcf. [12], the theoretical treatment is more involved. Moreover, 
the contribution of linearly polarized gluons is suppressed in %ci or Xbi production due to helicity conservation. An 
even stronger suppression arises for Xc2 and \b2 states, but as these can be produced from two on-shell gluons, they 
are theoretically cleaner and can help to cancel out uncertainties as we will discuss. 

Recently it was pointed out [13, 14] that in the production of Higgs bosons at LHC, there is in principle also a 
contribution from linearly polarized gluons inside unpolarized hadrons [15]. This affects the transverse momentum 
distribution in such a characteristic way that it provides a tool to determine whether the Higgs is a scalar or pseu- 
doscalar boson. As the distribution of linearly polarized gluons inside unpolarized protons is currently unknown (only 
a theoretical upper bound has been derived [15, 16]), no model independent predictions for the modification of the 
transverse momentum distribution could be given. For this purpose an independent experimental determination is 
required and several suggestions have been put forward, such as through azimuthal asymmetries in heavy quark pair 
production in deep inelastic electron-proton scattering (with the heavy quarks having large relative and absolute 
transverse momenta) [16] and photon pair production in proton-proton collisions [17]. Here we will discuss another 
promising method to access and extract the distribution of linearly polarized gluons inside unpolarized hadrons, 
namely that of (pseudo-)scalar charmonium or bottomonium production, which is more closely related to Higgs pro- 
duction and has the additional advantage of providing a way to map out the scale dependence. The latter is dictated 
by the framework of transverse momentum dependent (TMD) factorization, which for Higgs production including 
linearly polarized gluons has been discussed in Rcf. [18]. In a collinear factorization approach, polarized gluons are 
generated perturbatively [19-21], while in the framework of TMD factorization polarized gluons are already present 
at tree level through a nonperturbative distribution, here denoted by h 1 9 . It corresponds to an interference between 
+1 and —1 helicity gluon states that would be suppressed without transverse momentum. As a consequence, h 1 9 
shows up in transverse momentum distributions, either of pairs of particles or of single particles. It can modulate 
the angular part of the transverse momentum distribution, but it can also generate a term in the cross section that 
is independent of the azimuthal angle. This happens when two linearly polarized gluons, one from each hadron, 
participate in the scattering and the final state particle is a scalar or pseudoscalar. In this paper we will study how 
the distribution of linearly polarized gluons may affect the angular-independent transverse momentum distributions of 
±+ quarkonia, employing TMD factorization in combination with the nonrelativistic color-singlet quarkonium model, 
which we will justify from the framework of nonrelativistic QCD (NRQCD). NRQCD in combination with TMD fac- 
torization, has been used in Ref. [22] to study single spin asymmetries in quarkonium production in lcpton-nucleon 
and nucleon-nucleon collisions. It was observed that the asymmetry is nonzero in ep^ collisions only in the color-octet 
model, whereas in pp^ collisions only in the color-singlet model. In Ref. [23] the asymmetry for J/ip production in 
ep^ collisions was studied using the color evaporation model. 

In the color-singlet model (CSM) [24-26], the quarkonium bound state is produced in a color-singlet state from a 
heavy quark and antiquark with small relative momenta. Although the CSM fails to describe the large transverse 
momentum spectra of C = — vector states J/ip, tp(2S), and T (for a discussion cf. Refs. [5, 27]), it can be justified from 
NRQCD for C = + states, like for 775 in Ref. [9] . In the framework of NRQCD [28-30] , hadroproduction of quarkonia is 
described in terms of a double power series expansion in the strong coupling constant a s and in the velocity parameter 
v, such that (Mqv 2 ) 2 <C (Mqv) 2 <C Mq, with Mq the heavy quark mass. The magnitude of the velocity is given 
by the self-consistency condition v ~ o. 8 (Mqv), yielding v 2 ~ 0.3 for charmonium and v 2 ~ 0.1 for bottomonium. 
For details we refer to Refs. [28-30]. According to NRQCD, a heavy quark-antiquark pair can be produced at short 
distances not just as a color-singlet, but also in a color-octet configuration, which subsequently evolves into a physical 
quarkonium state by radiating soft gluons. The hadronization of the pair is encoded in long-distance, universal, 
matrix elements that can be characterized according to the velocity expansion in v. These matrix elements are not 
calculable perturbatively and have to be determined by a fit to the data. Lattice QCD estimates are available only 
for annihilation matrix elements [31]. NRQCD has had success in explaining many experimental observations, but 
it is not able to reproduce in a consistent way all cross sections and polarization measurements for charmonia, see 
Refs. [1, 2, 5] and references therein. Nevertheless, for the low transverse momentum part of the spectrum of C = + 
states, we expect it to provide a reasonable description, such as for the Xci,2 states discussed in Ref. [27]. For C = + 
states an NRQCD analysis [28] shows that color-octet contributions are (at least) order v 2 suppressed with respect 
to the singlet contributions. Earlier, in a study [12] of hadroproduction of Xcj in the fc^-factorization approach, it 
was also noted that the color-octet part is strongly suppressed and the color-singlet part dominates, especially at low 
transverse momenta. Also, according to Refs. [1, 2], color-octet contributions can certainly be neglected for C = + 
bottomonium, in agreement with the analysis of Ref. [9]. At large transverse momentum [32, 33] and central exclusive 
production [34] complications may arise, but these are beyond the scope of this paper. At low transverse momenta we 
consider it justified to restrict to the color-singlet contributions for our particular analysis, which is primarily aimed 
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FIG. 1: Feynman diagrams for the leading order gluon-gluon fusion process in the color-singlet model. 

at pointing out how the linear polarization of gluons enters in the production cross sections of C = + quarkonium 
states. 

II. OUTLINE OF THE CALCULATIONS 

We consider the process 

h(P A )+h(P B ) -+ QQ[ 2S+1 L { p]{q)+X , (1) 

where the four-momenta of the particles are given within round brackets, and we assume that the colorless heavy 
quark-antiquark pair is in a bound state described by a nonrelativistic wave function with spin S, orbital angular 
momentum L and total angular momentum J. The S, L, and J quantum numbers are indicated within the square 
brackets in spectroscopic notation, while the color assignments of the pair are generally specified by singlet or octet 
superscripts, (1) or (8). The squared invariant mass of the QQ state is M 2 = q 2 and M — 2 Mq up to small relativistic 
corrections. According to the CSM [24-26], we assume that the two quarks are produced in a hard partonic scattering 
with the same quantum numbers as the mesons into which they nonperturbatively evolve. Therefore, to lowest order 
in perturbative QCD, <D{a 2 ), one only has the gluon fusion process described by the two Feynman diagrams depicted 
in Fig. 1. This production mechanism is only significant in the low q T region with a mean q T determined by the 
intrinsic transverse motion of the gluon constituents. In practice, a fixed order calculation is not expected to describe 
very accurately the small transverse momentum spectrum [9] , but corrections and resummation can be considered at 
a later stage when it becomes relevant. Moreover, we expect that in the comparison of different quarkonia the effect 
of higher order corrections will matter less. Another aspect to keep in mind is that at small transverse momentum, 
but large xp, higher Fock states are expected to become relevant, such as considered in forward Higgs production in 
Ref. [35]. 

Transverse momentum dependent gluon distributions inside unpolarizcd hadrons are defined via a matrix element 
of a correlator of the gluon field strengths F^(0) and i< w (£), evaluated at fixed light-front (LF) time £ + = £-n = 0, 
with n being a lightlike vector conjugate to the four-momentum P of the parent hadron. Expanding the gluon 
four- momentum as p = x P + p T + p~n, the correlator can be written as [15] 



= Jf^ S ^'^ ^ (P| T r [^(O)^(C)] \P)\ 



LF 



— [a*" f 9 (r n 2 ) - ( P ^ I a» v P 
2x\ 9t h ^ P -> I M 2 +9t 2M? 



ff( X ,pi) [^f+g^^ ) ht 9 (x,pi) I (2) 
h h > 



where p 2 — —p 2 , — — P^n v jP-n — n^P v / P-n and Mh is the hadron mass. The unpolarized and linearly 
polarized gluon distributions are denoted by ff(x,p 2 ) and h 1 9 (x,p 2 ), respectively. A Wilson line, which makes the 
correlator gauge invariant, has been omitted in Eq. (2). Like other TMDs, also h x 9 (x,p T ) receives contributions from 
initial and/or final state interactions that render the gauge link process dependent. Therefore, in principle hf 9 (x,p 2 ) 
can receive nonuniversal contributions, despite the fact that it is T even, and its extraction can be hampered in those 
processes for which factorization does not hold, e.g. dijet production in hadron-hadron collisions [36-39]. For the 
present process we expect no problems with factorization since the C-even quarkonium bound state forms a color 
singlet. In this case there is only one color trace present and gauge links can always be factorized. 

Along the lines of Refs. [13, 39, 40], we assume that at sufficiently high energies TMD factorization holds and the 
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cross section for the process in Eq. (1) is given by 

dcr = y s ^3 2q0 j dx a dx b d 2 p aT d 2 p bT (2n) 4 5 4 (p a +p b -q) 



xTr I $ g (x a ,p aT )<l> g (x b ,p bT ) \-A (.9.9 ->■ QQ[ 2S+1 L ( ?]) ( Pa ,p b ;q)\ \ , (3) 

lv colors ) 

where s = (Pa + Pb) 2 is the total energy squared in the hadronic center-of-mass frame. Here and in following 
expressions, p a and p b are understood to have Pa = P b = 0- The amplitude A can be written in the form [24-26] 

A(gg^QQ[ 2S+1 L { ?]) (p a ,p b ;q) = J ^ Tr [0{ Pa ,p b] k) <j> LS -jj z {q, *0] , (4) 

with 2k being the relative momentum of the heavy quarks, see Fig. 1. In Eq. (4) (f>LS;JJ z (q, fc) is the Bcthc-Salpeter 
wave function of the produced bound state and 0(p a ,p b ; k) is calculated from the Feynman diagrams in Fig. 1 without 
including the polarization vectors of the initial gluons and the heavy quark legs, since they are absorbed into the 
definitions of the gluon correlators and of the bound state wave function, respectively. Therefore we have 



0^(p a ,p b ;k) = 2 5s 2 ^( 3 z;3j|l) \(T a T%rj^ 
id I [Pa 



j) b -j> a + 2$ + 2Mq „ 
7 



p b - 2kf - AM 2 



, (TbTa) . . , ta-t b + 2# + 2M Q 

+ (TT)y7 ( Pa -p b + 2k) 2 -AM 2 / /' (5) 

where the sum is taken over the colors of the outgoing quark and antiquark and T a are the SU (3) matrices normalized 
according to Tr(T Q T b ) = Ti(T b T a ) = 5 ab /2. The 577(3) Clebsch-Gordan coefficients, 

5 ij 

(3z;3j|l> = — = , (6) 



with N c being the number of colors, project out the color-singlet configuration. In the rest frame of the bound 
state, the relative momentum of the two quarks is small compared to their mass Mq, which justifies a nonrelativistic 
approach. Hence (pLS;.J.J z can be written as 

<pLS-,jj z (q,k) = 2ns(k°-^-\ Y,^LL z {k){LL z] SS z \JJ z )Vss z {d,k), (7) 

^ ' L Z ,S Z 

where ^LL z (k) is the eigenfunction of the orbital angular momentum L and the brackets denote the appropriate 
Clebsch-Gordan coefficients. The spin projection operator Vss z is defined in terms of the heavy quark and antiquark 
spinors, u (q/2 + k, s) and v (q/2 — fc, s), respectively, with s and s being their spins. More explicitly, 

,1 1 



VssMk) = ^2 S; 2~ slSSz)u (l + k ' S )^(l' k ' S ) 



1 



= 4M372 (t + W + M ) U ss z (-4 + n + M) + 0(k 2 ) , (8) 

with Hss z = 7 5 for singlet (S = 0) states and Uss z = $S z (q) for triplet (S = 1) states, where £s z (q) is the spin 
polarization vector of the QQ system. Substituting Eq. (7) into Eq. (4), one gets 

/j3i 
— ^^ LL Xk){LL z -SS z \JJ z )^[0{p a ,p b -,k)4> LS ,jj z {qM , (9) 



L Z ,S Z 



which can be expanded in powers of |fc| around k = 0. At this point it is convenient to separate the Fourier transform 
of the wave function *f>LL z (k) into its radial and angular pieces, 



/ 



^ e* kr * LLz (k) = * LLz (r) = R L {\r\)Y LLz (9, <p) , (10) 
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whcrer = (\r\,9,ip) in spherical coordinates, Rl(\t\) is the radial wave function and Y LLz (9, ip) is a spherical harmonic. 
This implies, in particular, 



[ ^ T 



i 



i?o(0) 



(11) 



Therefore, for S waves (L = 0, J = 0,1), the first term of the Taylor expansion of the amplitude is obtained by 
calculating the integrand in Eq. (9) at k — and, since this does not depend on k anymore, it can be factored out 
and one is left with the integral in Eq. (11). Making use of Eq. (8) as well, one gets 

A[ 2J+1 sf]{ Pa , Pb -q) = -L R (0)Tr[O( Pa ,p b ;0)Vss z (Q,Q)] > (12) 



(13) 
(14) 



4ir 

which, adopting the definition O(0) = 0{p ai Pb\ 0) with p b = q — p a , leads to the final results [26]: 



AfS^Kp^q) 



1 



4vVM 
1 



i? (0)Tr[O(0) ii + M)^} 
i? (0)Tr[O(0) (4 + M)ts,] 



4V7rM 

For P waves (L = 1, J = 0, 1, 2) i?i(0) = 0, so one needs to consider the linear term in k a in the expansion of Eq. 
(9). Again from Eq. (10) it can be shown that 



(15) 



where e£ {q) is a polarization vector referring to an L = 1 bound state, and R[(0) is the derivative of the P-wave 
(radial) wave function evaluated at the origin. One gets [26] 



A[ 2S+1 P^}( Pa ,q) = -i^I-B^O) J2 (lL z ;SS z \JJ z )el z (q) ±-Tr[0{p a ,q-k)V S sA<lM 

L z ,S Z 



(16) 



fc=0 



The numerical values of the radial wave functions and their derivatives at the origin are commonly obtained from 
potential models, or from leptonic decay matrix elements [41, 42]. For the singlet state, Eq. (16) reduces to 



A[ 1 P^}(p a ,q) = -i 



AttM 



#1(0) Tr 



O(0)^+elO a (0) 2 



M 



where we have defined 



(17) 



(18) 



In order to write the triplet amplitudes in an analogous form, we utilize the following relations for the Clcbsch-Gordan 
coefficients and the various polarization vectors [24, 25] : 



J2 (lL z ;lS z \00)e a Sz (q)4M) = \f\ (f P ~ ^f) > 

L Z ,S Z 

J2 (lL z ;lS z \U z )e a Sz (q)el(q) = -j- J\e, vpr7 g»« ? e\{q) 

L Z ,S Z 

]T (lL z ;lS z \2J z )e a Sz (q)el(q) = ef{q) . 



L Z ,S Z 



Hence, from Eq. (16), it is straightforward to express the triplet amplitudes as follows [24, 25]: 

i + M~ 



A[ 3 P^](p a ,q) 
A[ 3 P^]( Pa ,q) 



A[ 3 P^]( Pai q) 



AttM 



J?i(0)Tr 



:u; ( ()) + ( 7a 6 a (0)-^d«(0) 



^^(0)w4^(,)Tr 



7^(0)^ + 0(0)1^' 



AttM 



R[(0)e^(q)Tr 



7a 0/3(0) 



i + M 



(19) 

(20) 
(21) 
(22) 
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Here e v j is the polarization vector for a bound state with J = 1 and it obeys the usual relations (which hold also for 

e v Sz and e£j, 



el(q)q a =0, £ e% (q) e*£ (g) = -g a ^ + 



q a q si 

~m 2 



Q 



af) 



(23) 



while e a f is the polarization tensor for a spin 2 system and satisfies 

ef z (?) = 4» , e£ a («) - , q a ef{q) = , 
£j"((l)£jf(q) = \ [Q^Q^ + Q^Q va ] - \ Q^ v Q afj . 



(24) 



By substituting the explicit expressions of the operators O(0) and O Q (0), 



(25) 



derived from Eq. (5) into Eqs. (13), (14), (17) and (20)-(22), it is found that the only nonzero amplitudes are the ones 
corresponding to the 1 Sq (t]q) and 3 Pq,2 (xqo,2) states, where Q — c,b: 



A^S^p^q) 



-2i- 



g ab 



<N C VttM 5 



Ro(0) ^Papq., 



A^[ 3 P^](p a ,q) = -2i- 



tab 



9, 



A^[ 3 P^]( Pa ,q) 



-2i- 



R' (0) 



N r V 7rM 3 



M 



2 g^PapPaa ~ 9^9* ~ 9p9% 



(26) 
(27) 
(28) 



After squaring and averaging over the colors of the incoming gluons, the corresponding cross sections are calculated 
using Eq. (3). 

Here the cross sections will be written in terms of operator matrix elements appearing in NRQCD. The long-distance 
matrix element of relevance for the t]q is related to the t\q radial wave function at the origin by the following relation 

[28]: 



(0|O^( 1 5 )|0) = ^|ii (0)| 2 [l + O( t ; 4 )], 



(29) 



where v is the velocity of the heavy quarks inside the QQ bound state. The other color singlet matrix element, 
(0\Vi Q ( 1 S'o)|0), starts out order v 2 suppressed and will be dropped. The color-octet contributions 1 Sq 8 \ 3 s[ S \ 1 P^ 
are, according to the NRQCD counting rules, suppressed by v 4 ,v 3 , and v 4 , respectively [27, 28] with respect to the 
leading color-singlet term and will also not be included. 

Similarly, for the \QJ states the matrix element of relevance is [28] 



(O\0* QJ ( 6 Pj)\O) 



3N C 
2n 



|iii(0)| 2 [l + (V)], J = 0, 1, 2 . 



(30) 



This matrix clement is order v 2 suppressed with respect to the matrix element in Eq. (29) for t\q [29], but is the 
leading one for \QJ states. However, according to the NRQCD counting rules, there is also a color-octet contribution 
of the same order in v 2 . In Ref. [27] it is shown that its matrix element is nevertheless suppressed, as reflected by 
the ratio R x ® = M^(0$ QJ ( 3 Si)) / (0* QJ ( 3 Pj)) w v°/(2N c ) w 0.17, where now the suppression comes from the color 
factor N c . The suppression of the color-octet contribution in \c production was also noted in Ref. [12]. 

Expressing the CSM results in terms of the leading color-singlet NRQCD matrix elements in Eqs. (29) and (30), 
yields the production cross sections presented in the next section for J = and J = 2 states. For J = 1 there is no 
contribution. 
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III. QUARKONIUM PRODUCTION CROSS SECTIONS AT SMALL TRANSVERSE MOMENTUM 



The calculation as explained in the previous section yields our main results for the production cross sections of 
J = and J = 2 C = + quarkonia: 



2ttV 



dyd 2 q T 9M 3 s 



(o\orCs )\o)c[ffff] [l-R(q 2 T )] , 



(31) 



dg(xQo) = 8 ir 3 a 2 
dyd 2 q T ~~~ 3 M 5 s 



(0\O* Qa ( 3 P Q )\Q)C{f?f?} [l + R(q 2 T )] 



(32) 



da(xQ 2 ) = 32 Tr 3 a 2 s 
dyd 2 q T 9 M 5 s 



(o\o^ 2 Cp 2 )\o)c[fffp 



(33) 



with N c — 3 and corrections of 0(q 2 /M 2 ) and C(u 2 ). Here y is the rapidity of the produced bound state along the 
direction of the incoming hadrons and R{q%) is defined as the ratio 



R(q 2 T ) 



c 



C [/?/?] ' 
The TMD convolutions are given by 

C[ W ff] EE 



d p aT / dp bT S (p aT +p bT - q T )w(p aT ,p bT ) f(x a ,p aT ) f(x b ,P bT ) 



with the transverse momentum weight 



1 



w 



2M h 



(PaT-Pbr) 2 - ^pl T pl T 



and, neglecting again terms of O (q 2 /M 2 ), 

M „ 



_ M_ y 



x b = —= e 



(34) 



(35) 



(36) 



(37) 



The results in Eqs. (31)-(33), integrated over q T , are in agreement with the ones published in Ref. [26] and ob- 
tained without the inclusion of intrinsic transverse momenta of the incoming gluons. Furthermore, we note that 
/ d 2 q T (q 2 ) a C[wH h^ 9 h^ 9 ] = for a = 0, 1, model independently. The case of a — implies that linearly polarized 
gluons do not affect the q T -integrated cross section. This means that if one normalizes the differential cross sections 
to the integrated ones, the uncertainty from the hadronic matrix elements (0|Ci |0) cancels out. The case of a = 1 
implies that R must have at least two nodes in q T , i.e. it has to flip sign at least twice, with 1 i?(qj = 0) ^ 0. This 
distinctive double node feature can be used to experimentally demonstrate the presence of the linear polarization of 
the gluons. 

It should be noted that the sign difference in the R{q%) term in Eqs. (31) and (32), for different parities of the 
quarkonia, is completely analogous to the case of pseudoscalar and scalar Higgs boson production, respectively, 
considered in Ref. [13]. There is a clear difference between the quarkonium and Higgs production cases though. The 
existence of a scalar Higgs boson does not automatically imply the existence of a pseudoscalar Higgs boson (e.g. from 
a supersymmetric extension of the Standard Model 2 ), and even if it does exist, it is far from given that it couples 
to Standard Model particles in such a way that this sign change could be observable. In contrast, both types of 
charmonium and bottomonium states follow from the same theory and have been established experimentally, thus 
offering a definite possibility to cross check the predicted features, i.e. the double node modulation of the cross section 
and the sign difference. 



1 If 9 is of definite sign and thus has no nodes itself, then R(q%. = 0) > 0. 

2 The recently observed Higgs boson signal at ATLAS and CMS cannot be accounted for by the CP-odd Higgs boson of the MSSM or 
NMSSM [43]. 




p T (GeV) p T (GeV) 

FIG. 2: Gaussian distributions for ff and hf 9 [divided by ff (x) ] in units of GeV~ 2 as functions of p T for two different values 
of (p%) with r — 2/3 (left) and r — 1/3 (right). 



The XQ2 cross section does not receive a contribution from linearly polarized gluons, as it would require a helicity 
flip of 4 units, which is heavily suppressed. As a consequence, to leading order in v 2 the ratio of xqo to \Q2 cross 
section yields a direct probe of the ratio R(q 2 ), as it simply equals |(1 + R(q 2 )). Both the uncertainty from the 
hadronic matrix elements (0|Oi|0) [as follows from Eq. (30)] and from the unpolarized gluon TMD cancel out in this 
ratio. For the Higgs boson case, no analogous ratio can be formed. Another difference that may turn out to be crucial 
is due to the scale dependence. The contribution from linearly polarized gluons is expected to decrease with increasing 
scales, just as other TMD effects, such as the Sivers single spin asymmetry [44-46]. Therefore, for quarkonia with 
mass of order 3 or 10 GeV, the effect of linearly polarized gluons is expected to be larger than at a Higgs mass scale 
of 125 GeV. Of course, if both quarkonium and Higgs production can be studied, then that would allow for a valuable 
check of the scale dependence expected from TMD factorization. 

It may be interesting to compare these findings to the expectations for the production cross sections in the color 
evaporation model (CEM) 3 . In the CEM the probability to form a physical (colorless) quarkonium state is assumed to 

be independent of the color of the QQ pair, which means that gg QQ^^L^p] is assumed to produce a colorless 
quarkonium state by means of soft gluon emissions that do not affect the momenta or angular momenta. The cross 
sections for gg — > QQ^S^], gg — > QQ^Pq ], and gg — > QQ^P^], then follow the same pattern as in Eqs. (31)-(33), 
namely R enters with a positive (negative) sign for a P-even (odd) bound state, while R = when J — 2. When 
integrated over q T the results for these octet contributions, agree with the ones published in Ref. [30]. Therefore, we 
expect inclusion of color-octet contributions not to affect our conclusions: they are suppressed in the NRQCD based 
calculation and in the CEM they yield the same results as in the CSM. In the CEM sometimes the probability is also 
assumed to be independent of the spin of the QQ pair [47] . In that case no difference between the various XQJ states 
should exist. If R{q 2 ) is found to be nonzero in XQO, this version of the CEM would be ruled out. 

If the initial gg pair is in a color singlet, the resulting QQ pair is necessarily in a C = + state [48], but for the 
color-octet case it can also yield C — — states. This means that in the CEM also C = — quarkonium states can be 
produced in the 2 — » 1 scattering process gg — > QQ[ 2S+1 L^j ']. However, for J = 1, such as for J/ip or T production, 
there is no contribution from gg — > QQpS^ 8 '] because of the Landau- Yang theorem, like for XQi- Therefore, also 
from CEM considerations we conclude that J/ip, tp(2S) or T production will not be useful for studies of hf 9 at small 
transverse momentum. 

At order a 3 1 ±:!= states can be produced through the subprocess gg QQg, but the contribution from hf 9 to 
the cross section is found to be power suppressed in these and all other quarkonium cases. Therefore, we do not 
expect any unsuppressed contribution from hf 9 in C = ± quarkonium production at large transverse momentum in 



3 For a detailed comparison between the CEM and NRQCD cf. [27]. 
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hadronic collisions. Of course it may be possible that for example gg — > Xco J/ip 7 produces a contribution for 
J ftp's with larger transverse momentum, but this will only be significant on the intermediate J — resonance which 
itself has small transverse momentum. Such cases simply follow the pattern of the results obtained in this paper. 
Note that in this example the final state photon (or gluon jet) momentum needs to be measured explicitly, otherwise 
no unsuppressed contribution from h t 9 arises. 



IV. NUMERICAL RESULTS 

In this section we estimate the possible size of the contribution from h^ 9 to the quarkonium production cross 
section, calculated at order a 2 within the framework of the color-singlet model. The effect will be maximal when the 
unknown function h^ 9 saturates its model-independent positivity bound [15], 

P7 ' \hi 9 (x,p 2 T )\ <f!(x )P 2 T ), (38) 



2MZ 



which is valid for all values of x and p T . Adopting a standard approach, we assume that the gluon TMDs have a 
simple Gaussian dependence on transverse momentum [13], 

ff( x , P l) = lM exp (ri\ , (39) 



""(Pt) V (P 

where f 9 (x) is the collinear gluon distribution and the width (pj) will in general depend on the energy scale, which 
is set by the quarkonium mass M. Moreover is taken to be independent of x. The bound in Eq. (38) is satisfied 
(but not everywhere saturated) by the form 

L -U, 2^ M 2 f?(x)2(l-r) / 1 pi \ 

hf 9 (x,p 2 T ) = \)J exp 1 - - ff- . (40) 



Tr(p 2 ) 2 r \ r (p, 

In Fig. 2 the p T dependence of ff and h^ 9 is shown for two values of the input parameter r (r = 2/3 and r = 1/3) and 
two values of the Gaussian width: (p^,) = 0.25 GeV 2 and (p 2 ^) = 1 GeV 2 , which are considered relevant for charmonium 
and bottomonium, respectively. Recall that the expressions given in the previous section apply for q 2 <§; M 2 . This 
means that for q 2 ~ M 2 we expect the approach and in particular the Gaussian falloff to be unrealistic, rather one 
expects a power-law falloff of hf 9 that would result from a perturbative collinear treatment [18]. In practice, this 
difference may not matter much, as the cross section and the modulation from linearly polarized gluons are expected 
to be quite small already in this tail region. 
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FIG. 4: Transverse momentum distributions for the different quarkonia, using the distribution functions in Eqs. (39) and (40) 
with r — 2/3, for two different values of (p%). 



The convolution integrals that appear in Eqs. (31)-(33) can be evaluated analytically: 
f?(x a )f?(x b ) _ ( 1 q 2 



9 m _ 



C[whf 9 hf 9 

Their ratio, 



-Lfl 1,-1-31 _ fl( X a)fl( x b) 



exp 



2 



4ir(p 2 



r(l - r) 2 



1 9? i (g 2 ) 2 



r (p 2 ) 8r 2 (p 2 ) 2 



exp 2 



1 9? 



2r <p2) 



1 



R(r,q 2 T ) = -r(l-r) 2 



1 g 2 1 {q 2 T f 



r (p%) 8r 2 (p 2 )^ 



exp 2 



1-r 9 2 



2r (pi) J ' 



(41) 
(42) 

(43) 



measures the relative size of the contribution by linearly polarized gluons. The behavior of R as a function of q T , 
presented in Fig. 3 for r = 2/3 (left panel) and for r = 1/3 (right panel), is characterized by the expected double node 
structure of C[w h\ 9 h\ 9 \ [13]. 
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FIG. 5: Same as in Fig. 4, but with r = 1/3. 
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The transverse momentum distributions of t\q and Xqo,2 are therefore given by 



1 dcr(xQo) | 7-)/ „2M 1 



o{XQ) Aql 2(p2) -V 2 (p2) 



,2 



[* + o73v«p --^ , (45) 



1 dtr(xQ2) 



.2 



exp -TT^T - (46) 



and they are shown in Figs. 4 and 5 for different values of the input parameter r and the Gaussian width (p%). As 
mentioned, the distributions for r] c ^ and Xc,bo are similar to the ones for a pseudoscalar and a scalar Higgs boson, cf. 
[14]. Unfortunately, there are no data available yet to fit or bound the modulation, neither for charmonium nor for 
bottomonium. Hopefully, LHCb will be able to provide such data in the near future. Existing data for Xc production 
([49] and more recently, from Tevatron [50] and LHC [51]) is from Xc —> J/i> 1, which is predominantly Xci and Xc2- 
Other decay channels, with a larger Xco component, should be considered instead. 

To obtain numerical estimates for the production cross sections, we have employed a Gaussian model. Since the x 
values of the gluons in one of the two protons can be rather small at LHC, one may also consider to calculate the cross 
sections within a small-x formalism. It has been shown that the distribution of linearly polarized gluons in various 
small-x treatments, the color dipole model [52, 53] and a color glass condensate model [52, 54], is quite sizable, in fact, 
can even saturate its bound. Therefore, we do not expect small- x effects to necessarily lead to smaller modulations 
compared to those obtained here for the Gaussian model. 

V. SUMMARY AND CONCLUSIONS 

In this paper C = + quarkonium production in hadronic collisions has been studied within the framework of 
TMD factorization in combination with the nonrelativistic QCD based color-singlet quarkonium model. It is shown 
that the transverse momentum spectra of scalar (XcO,Xbo) and pseudoscalar quarkonium (rj c ,rjb) production arc 
modified in different ways by linearly polarized gluons that can be present inside unpolarized hadrons. The scalar 
spectra receive a modulation R(q%,), whereas for pseudoscalars it is of opposite sign, — R(q^,). The modulation itself 
displays a characteristic double node structure, the magnitude of which is not yet known however. The advantage of 
using (pseudo)scalar final states is that no angular analysis needs to be performed to probe the effects of the linear 
polarization of gluons. Their effects on the production of higher angular momentum quarkonium states are strongly 
suppressed. As a consequence, comparison of xqo and XQi production can help to cancel out uncertainties. Similar 
results are not obtained for C = — states, as they turn out not to be sensitive to the linear polarization at leading 
order. We thus conclude that C = + charmonium and bottomonium states offer an excellent way to access and extract 
the distribution of linearly polarized gluons inside unpolarized hadrons and to check the predicted features: the double 
node structure and the sign difference between scalar and pseudoscalar production. Together with the analogous study 
in Higgs production at LHC, quarkonium production can moreover be used to test the scale dependence of the linearly 
polarized gluon distribution over a large energy range. Experimental opportunities to study these new proposals are 
currently offered by LHCb and the proposed fixed-target experiment AFTER at LHC. 
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